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The decay of massive neutrinos to final states containing 
only invisible particles is poorly constrained experimentally. 
In this letter we describe the constraints that can be put on 
neutrino mass and lifetime using CMBR measurements. We 
find that very tight lifetime limits on neutrinos in the mass 
range 10 eV - 100 keV can be derived using CMBR data from 
upcoming satellite measurements. 
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If the mass of any given neutrino species is larger than 
roughly 100 eV it is necessarily unstable with a lifetime 
shorter than the Hubble expansion timescale JIJ] . A mas- 
sive neutrino could in principle have many different al- 
lowed decay modes, but if the final state contains elec- 
tromagnetically interacting particles the decay is ruled 
out by observations [||, unless the lifetime is very long 
as in the scenario proposed by Sciama ||. However, if 
the final state contains only "invisible" particles (mass- 
less neutrino species, majorons, etc.) it is very difficult 
to constrain the decay by looking for the decay products, 
and in this case it becomes important to use indirect ar- 
guments. In this letter we shall look at this latter type 
of decay, where the decay products cannot be directly 
observed. 

By using arguments pertaining to big bang nucleosyn- 
thesis one can rule out a large chunk of the mass-lifetime 
parameter space j|. In essence, the argument is that if 
the mass is too large and the lifetime too long, there is 
too much energy density present in the universe during 
nucleosynthesis which in turn leads to an overproduction 
of helium Q . If the heavy neutrino is meta stable during 
nucleosynthesis, one can rule out a neutrino mass larger 
than a fraction of an MeV, but if the neutrino is unstable 
on the timescale of nucleosynthesis (r < 10 3 s), the argu- 
ment becomes more complicated, and certain lifetimes 
are allowed even for heavy neutrinos. Altogether, nucle- 
osynthesis can be used to constrain lifetimes for masses 
larger than roughly 0.1 MeV 

In the present letter we focus on another method that 
can be used for the same purpose. The cosmic microwave 
background radiation (CMBR) is created at the last scat- 
tering surface for photons. The CMBR is anisotropic at 
the 10 -5 level, a fact long predicted but only detected in 
1992 by the COBE satellite [g. The physical conditions 



prevailing at the epoch of last scattering are imprinted in 
these anisotropies, meaning that one can in principle de- 
termine the physical content of the universe at photon de- 
coupling by measuring the anisotropy spectrum [Q . Pho- 
tons decouple from the cosmic plasma at a temperature 
of approximately 1 eV and if conditions are changed (for 
example by a heavy, decaying neutrino) prior to this, it 
is possibly noticeable in the CMBR anisotropics. There- 
fore it should be possible to use such arguments to con- 
strain decays happening as late as at a temperature in 
the eV range. The influence of decaying neutrinos on 
structure formation have been considered many times in 
the literature H|| and it is also well known that decay- 
ing neutrinos may influence the CMBR power spectrum 
However, such calculations have not focussed on the 
limits that may be put on neutrino masses and lifetimes. 
Previously, the data have not been sufficient, but within 
the next few years we can expect to have extremely fine 
data for the power spectrum of microwave anisotropies, 
measured by the upcoming satellite missions MAP and 
PLANCK p0[ . This is the primary motivation for dis- 
cussing the neutrino mass and lifetime limits that can 
possibly be obtained from CMBR data. 

We shall here assume that the decaying neutrino is 
the r neutrino and that the two lighter neutrinos are 
massless (that is, they have masses much smaller than 
the photon decoupling temperature). The same type of 
argument could also be used for an unstable muon neu- 
trino, but if the muon neutrino is massive and decays on 
this timescale, it would also imply a heavy and unstable 
tau neutrino unless there is some kind of inverted mass 
hierarchy. Having a scenario with two decaying neutrino 
species complicates the matter unnecessarily although 
such a scenario could have interesting consequences [[llj . 
The decay we study can therefore be written generically 
as 
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where D indicates the final state (yvv, vcj), etc.) config- 
uration. We only look at decays taking place after Big 
Bang nucleosynthesis; that is, we assume m Ur A T(r) < 
0.1 MeV. 

If the decay final state contains only massless particles 
and we ignore possible inverse decays as well as relativis- 
tic corrections the equations describing the decay sim- 
plify a lot. We follow for instance Dodelson, Gyuk and 
Turner |l2| in taking this approach. In that case, the 
cosmological evolution equations have the form 
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where To and po indicate the temperature and energy 
density of a standard massless neutrino species. The 
subscript D indicates the decay products and r is the 
lifetime of v T . We assume that the decay is finished well 
before the last scattering (in practice we assume that nei- 
ther the heavy neutrino mass nor the temperature corre- 
sponding to r are smaller than 10 eV). This means that 
the decay only perturbs the CMBR through the change 
in relativistic energy density at last scattering. By using 
these assumptions the decay simplifies further because 
the final relativistic energy density depends on m Ur and 
t only through the combination r. One can define a 
decay "relativity" parameter, a, as [jy] 
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the decay being relativistic roughly for a ^ 1 and non- 
relativistic for a > 1. We follow the standard practice 
in expressing the relativistic energy density in equivalent 
number of massless neutrino species ]13| 



N v = (2p + p D )/ Po . 



(4) 



Within our approximation, N v — > 3 for relativistic de- 
cays, since in that case all that really happens is a 
reshuffling of relativistic energy density between different 
species. If one includes inverse decays, relativistic decays 
proceed in equilibrium so that the decay does not fin- 
ish before the heavy neutrino goes non-relativistic. This 
means that even for strongly relativistic decays there will 
be some contribution to the energy density of the decay 
products from the rest mass term, so that N v > 3 in this 
case. The decay formalism is much more complicated, 
however, since one then has to use a specific final state 
(depending on the number of final state particles and 
their quantum statistics) and solve the full Boltzmann 
equation For this reason we just use the simplistic 
approach because the main point of the letter is not to de- 
rive precision limits but rather to point out the strength 
of this type of argument. 

For very non-relativistic decays one can also arrive at 
an analytic approximation for N v . Assuming that the 
cosmic energy density just prior to decay is completely 



dominated by the heavy neutrino and that the decay pro- 
ceeds instantaneously at t = r, one gets N u ~ 0.52 a 2//3 . 
Altogether we can approximate N v for any value of a 
with the function 
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0.52 a 2/3 , 



(5) 



which has the correct asymptotic behaviour. In Fig. 1 
we have plotted both N v from the numerical solution 
of Eq. (Q) and from Eq. (J|). In the remainder of the 
letter we shall use our analytic approximation, Eq. (|^), 
to simplify calculations. 
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FIG. 1. The equivalent number of neutrino species pro- 
duced by decay. Both the solution obtained by use of Eq. (^) 
and our analytic approximation, Eq. (^), are shown. 

Note here, that we have only included the standard 
three neutrino species and the decay products in this cal- 
culation. In principle there could be energy density con- 
tributions from other more exotic particles, which would 
add some constant value to N v . We shall not deal fur- 
ther with this possibility here, only note that if this is 
the case, the allowed value of a corresponding to a given 
N v is shifted to lower values. 

The microwave background anisotropies are usually 
calculated in terms of C/ coefficients that relate to the 
fluctuation spectrum as C; = (|a/ m | 2 ), where a; m can 
be expressed in terms of the temperature fluctuations as 
T{8, 0) = J2im ai m Yi m (&, 0)- At present, CMBR data are 
quite insufficient to put any constraints on the equivalent 
number of neutrino species. The main reason is that the 
only real precision data at present are the COBE data 
||. This data only gives the power spectrum at I < 20, 
where the spectrum is flat and featureless. However, as 
previously mentioned this can be expected to change dra- 
matically within the next few years, where two different 
satellite missions are expected to measure the CMBR 
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power spectrum to high precision out to values of I be- 
yond 1000 This will allow for determination of al- 
most all of the relevant cosmological parameters to great 
precision, including the number of light neutrino species 
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FIG. 2. CMBR power spectra for models with different val- 
ues of a. The spectra have been normalised to the quadrupole 
coefficient, 6C2. The underlying model is a standard CDM 
model with Q — 1, fij, = 0.05 and h — 0.5. 

Several different groups have performed calculations of 
the precision to which one can measure different param- 
eters, depending on the quality of the data [Jl4],[l5|. It 
turns out that the CMBR power spectrum is very sen- 
sitive to N v because varying the relativistic energy den- 
sity also means shifting matter-radiation equality. In the 
matter dominated epoch photons see a constant gravita- 
tional potential (in time) after they are emitted from the 
last scattering surface, at least in the linear approxima- 
tion. If the universe is not completely matter dominated 
at last scattering the potential is not constant and this 
in turn increases the CMBR anisotropy, an effect known 
as the early ISW effect 0. If the amount of relativis- 
tic energy density is increased, last scattering will occur 
closer to the radiation dominated epoch, meaning that 
especially the first Doppler peak is increased in height. 
In Fig. 2 we show the CMBR power spectrum for differ- 
ent values of a. Clearly, the first Doppler peak is very 
sensitive to changes in a, exactly as expected. This high 
sensitivity to N v means that N v can be measured fairly 
precisely by CMBR experiments. As shown in Ref. [ jL4j , 
experiments like MAP and PLANCK should realistically 
be able to measure AN V to within 0.3 without any prior 
knowledge of other cosmological parameters, and to much 
greater precision assuming some prior knowledge of other 
parameters. Note here that by using nucleosynthesis ar- 
guments one can also constrain N v to be below roughly 



4, but this argument only holds for temperatures above 
the range we are interested in. In our case we assume 
that the decay takes place after nucleosynthesis so that 
during nucleosynthesis we have N v = 3. 

Now, if we assume a given mass for i/ T , we can relate 
AN V to an allowed lifetime interval At for any value of 
N v . Doing this gives 

± { N„ ± AN„ - 3 \ 3/2 _ 2 

r ( y ) I — rig — J 171 ^ ( ' 

where t + indicates the maximum and t~ the minimum 
allowed lifetime. 
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FIG. 3. Allowed lifetime intervals for neutrino decay, plot- 
ted for different values of measured N v . The three intervals 
correspond to N u = 3.5, 5 and 10; the interval with the longer 
lifetime corresponding to the higher N v . 

In Fig. 3 we have shown the allowed lifetime inter- 
vals for a decaying tau neutrino as a function of neutrino 
mass for different measured values of N v and assuming 
different precision. Especially for non-relativistic decays 
one can constrain the lifetime to lie within a very narrow 
band. Note that measuring N v ^ 3 is of course not a 
confirmation of neutrino decay, although that would be 
a possible explanation for such an effect. As such, the 
CMBR arguments (like nucleosynthesis arguments) can 
only serve to rule out certain regions of parameter space, 
not provide an unambiguous detection. 

In conclusion, we have shown that the CMBR is a very 
sensitive tool for constraining neutrino decays, and with 
new CMBR data it should become possible to narrow 
down the allowed lifetime intervals for neutrino masses 
in the range of 10 eV to 100 keV. The constraints coming 
from this type of argument are applicable to any type of 
neutrino decays to massless, non-interacting final states. 
What is constrained is really the amount of relativistic 
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energy density at last scattering, and in that regard this 
type of argument is very similar to nucleosynthesis con- 
straints. What is different with the CMBR argument is 
that it pertains to a completely different mass-lifetime 
region. Here, we constrain masses in the region below 
0.1 MeV, so that by putting together CMBR and nu- 
cleosynthesis arguments essentially all masses down to 
O(10eV) are covered, since nucleosynthesis already deals 
with more or less all masses up to the maximum exper- 



imentally allowed mass of 24 MeV |16|. Thus, although 
CMBR data cannot at present constrain neutrino decays 
to invisible final states it is to be expected that in the very 
near future this will change. One should then be able to 
constrain such neutrino decays very strongly with new 
data. 
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